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Abstract 

Information theoretic geometry near critical points in classical and quan¬ 
tum systems is well understood for exactly solvable systems. Here we show 
that renormalization group flow equations can be used to construct the 
information metric and its associated quantities near criticality, for both 
classical and quantum systems, in an universal manner. We study this met¬ 
ric in various cases and establish its scaling properties in several generic 
examples. Scaling relations on the parameter manifold involving scalar 
quantities are studied, and scaling exponents are identified. The meaning 
of the scalar curvature and the invariant geodesic distance in information 
geometry is established and substantiated from a renormalization group 
perspective. 
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1 Introduction 


Information geometry provides a unique arena where geometric notions are ap¬ 
plied to physical systems, often leading to new and alternative insights into the 
physics of classical and quantum phase transitions. A Riemannian metric defined 
on the parameter space for classical systems or the space of coupling constants 
for quantum systems define a distance on the parameter manifold (PM) [T], [2j. 
Geometric properties of this distance translate into useful physical quantities to 
understand phase transitions. Although this method is primarily used to study 
second order continuous transitions, first order phase transitions can also be 
treated in the geometric framework 0.0- 

Geometric methods have often been applied to statistical systems that are 
solvable. Namely, one calculates the metric and its associated quantities ana¬ 
lytically using an equation of state or a solvable Hamiltonian, and studies their 
limiting behaviour as one approaches criticality. For example, in classical liquid- 
gas or magnetic systems, one can use the Van der Waals equation of state to 
analytically compute the metric, or use experimental data (based on multi pa¬ 
rameter fits to equations of state) to compute the same [3], [3]. In the context of 
quantum systems, one normally alludes to a complex quantum geometric tensor 
whose real (symmetric) part is the Riemannian metric on the parameter mani¬ 
fold and its imaginary part is the Berry curvature. In exactly solvable quantum 
systems, a knowledge of the ground state leads to the metric [5j. 

Scaling analyses of metrics in quantum phase transitions were first performed 
in the important work of [6]. Here, the authors provided an integral representation 
of the quantum geometric tensor in terms of imaginary time correlation functions 
and hence were able to extract information regarding the scaling properties of the 
metric tensor (see also [7] for related work regarding the scaling behaviour of the 
Gaussian curvature in the context of the XY spin chain model). The broad issue 
that we address in this work is whether there is a generic way of understanding 
the geometry of phase transitions, both classical and quantum, particularly in 
situations where an exact solution to the model (or an equation of state) may 
not be available. 

To this end, we first note that there are indications that as far as geometry 
is concerned, the descriptions of classical and quantum phase transitions might 
be very similar. Let us briefly elaborate on this by focusing on two dimensional 
PMs which will mainly be of our interest here. Important in the study of any 
geometric setup are scalar invariants in these. These invariants, which are coor¬ 
dinate independent, provide an invariant characterisation of a curved manifold. 
As is well known, in two dimensions, the curvature scalar (or the Ricci scalar) 
completely characterises the curvature. Associated to this is the scalar expansion 
parameter (to be elaborated towards the end of this section) which measures how 
geodesics (which are analogues of straight lines in curved spaces) converge (or 
diverge) towards a point in the PM [Bj. Further, the line element, identified with 
an affine parameter that measures infinitesimal distances along geodesics provide 
a third scalar quantity. As pointed out in [9], Relations between these scalar 
quantities reveal universal behaviour in classical and quantum phase transitions. 
Namely, the scaling behaviour of the Ricci scalar and the expansion parameter 
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with the affine parameter near criticality is universal for these, in any two dimen¬ 
sional PM, under the assumption that the scalar curvature diverges at criticality 
as a power law. This hint of universality naturally leads one to suspect that there 
might be a generic way to compute metrics on the PM, at least near criticality 
and we indicate how this can be achieved, by using ideas from scaling symmetries 
near critical points. 

As far as classical phase transitions are concerned, the usefulness of this 
method is that we are able to compute the information metric in varied sce¬ 
nario, which we believe were hitherto not known. For example, as we elaborate 
upon later, the metric for Ising type models close to, and at, four dimensions 
computed from our method shows interesting and non-trivial behaviour of scalar 
invariants in information geometry, consistent with the physics near their fixed 
points and we will show that the scaling behaviour of the Ricci scalar acquires 
logarithmic corrections in this example. 

Although most of this work deals with classical phase transitions, we apply 
our method to one example in the context of zero temperature quantum phase 
transitions and compute the metric using scaling arguments, and show that we get 
consistent results. Our method here should be contrasted with the one developed 
in [6]. As we have mentioned, the latter used the imaginary time correlation 
function to derive the scaling relations of the metric in the context of quantum 
phase transitions. Here, we will directly appeal to the renormalization group (RG) 
equations to construct the metric and show that this gives sensible results. In a 
sense this complements the approach of [6], while retaining universal applicability. 

It is useful to point out here the conventional definitions of the information 
metric in classical and quantum systems. In the former case, Ruppeiner’s defini¬ 
tion of the metric [T| reduces to derivatives of a thermodynamic potential. For 
example, if we consider the entropy density s, then the line element and the 
metric on the PM are defined via 



( 1 ) 


where x a , a = 1,2, denotes the internal energy and the particle number (both 
per unit volume) and the Boltzmann’s constant is set to unity. These are the 
co-ordinates on the parameter manifold in the “entropy representation.” The 
metric can also be computed in various other representations and a full list of 
such metrics is available in Table II of [I]. We mention here that an alternative 
definition of the metric had been advocated earlier in the seminal works of Wein- 
hold [10]. That the two definitions are related by a conformal transformation is 
well known. 

On the other hand, the information metric in quantum systems is defined 
by considering two infinitesimally separated quantum states (in the parameter 
space) and computing 


\ip (x + dx) — i/j (x) | 2 = (d a ip\db^)dx a dx b = a a bdx a dx b 


( 2 ) 


where x a (collectively denoted as x in the l.h.s of eq.(]2])) denotes the parameters 
on which the wave function i[} depends on, and d a is a derivative with respect to 
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x a . The approach of p] is to construct, from a a b (which are not gauge invariant), 
a gauge-invariant metric tensor given by 


9ab = Oiab ~ PaPb\ /3 a = (x) \d a ^ (x)) 


( 3 ) 


Here, g a b is the metric induced on the PM from the natural structure of the 
Hilbert space of quantum states. Eqs. (JU and (j3]) are the standard definitions 
of the classical and quantum mechanical Riemannian metrics on the parameter 
manifold. Our approach here would be to compute these metrics in the critical 
regime, without using Eqs. (Op) and (EJ) directly. 

For this purpose, we use an existing notion in the literature, namely the ge¬ 
ometric equivalent of scale invariance near a fixed point. To the best of our 
knowledge, such a proposal first appeared in m- In this work, the scaling re¬ 
lations for classical liquid-gas phase transitions are recovered from a geometric 
perspective. Related work has appeared in the literature in the context of quan¬ 
tum field theory [12], [13], [H], [TH] . [16]. [T7j and statistical mechanics [T8] . The 
main idea in these works is that the renormalization group flow equations deter¬ 
mine the so called homothetic vector fields, which are mathematically related to 
scale invariance. 

Specifically, if K° is a homothetic vector held on a manifold with metric g ab 
(we will mostly consider two dimensional manifolds so that a,b = 1, 2), it satisfies 
the condition Cxg a b = Dg a b, where C-K9ab is the Lie derivative (see, e.g section 
1.4 of [19]) of the metric along a curve whose tangent is K a . Also D is a constant 
that we will identify with the spatial dimension of the system (note that this is 
different from the dimensionality of the parameter manifold, which will mostly 
be two in this paper). This equation reduces by standard manipulations to the 
condition K a . b + K b]a = Dg a b where a semicolon denotes a covariant derivative 
(defined in the next section) on the parameter manifold (again, see section 1.4 
of [T9]). This last equation can be alternatively written in a simpler form as 


( 4 ) 


9ac^ t b T 9bc-k -T 9ab,cK Dg a b j 


where the comma indicates an ordinary derivative with respect to the coordi¬ 
nate label that follows it, and repeated indices imply a summation (which will 
always be the case in this paper). To fix ideas, let us consider a textbook ex¬ 
ample, the two dimensional Euclidean space. This is flat space, with coordinates 
(x,y), and metric g ab = diag(l,l). Considering the vector field Ii a = (x,y), 
which represent the tangent at any point on the flat manifold, it is seen that 
K a -b + Kb-,a = diag(2, 2), thus confirming Eq.(j4]). This is expected, since two di¬ 
mensional Euclidean space is flat, and looks the same at any length scale. Curved 
manifolds, which will be of our interest here are more challenging to deal with. In¬ 
deed, properties of homothetic vectors (when they exist) are of great importance 
in General Relativity and Cosmology. Here, we will apply this to information 
geometry of phase transitions. 

In this paper, we will use the notion of scale invariance of the parameter 
manifold near a critical point. Our starting point is a metric on the parameter 
manifold of a system, i.e the space of the coupling constants in the theory. This 
metric is assumed to be a priori unknown. Near criticality, following [XT] , [12], we 
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demand that the beta functions of the theory are the components of a tangent 
vector field which is homothetic. From Eq.flJJ), we then get a set of coupled 
partical differential equations for the components of the metric. These equations, 
if solvable, will lead to solutions of the metric on the PM, without a detailed 
knowledge of the full solution of the system. 

For classical systems where, in the conventional approach of Ruppeiner [I], 
the metric components are defined via derivatives of the free energy (or entropy) 
and are related to response functions, this has been demonstrated in mi. Here, 
up to linear order in RG, it was found that Eq.Q implied that the metric compo¬ 
nents are generalized homogeneous functions. Euler’s theorem was then invoked 
to read off the scaling behaviour of the metric components. In this paper, we 
will consider situations where this may not be possible, and solve for the metric 
components directly from Eq. fjTjl . Using scale invariance of the parameter man¬ 
ifold, our method should be viewed as a tool for obtaining the geometry of any 
system, sufficiently near to criticality. This non-trivially generalizes the analysis 
of m by providing an universal approach to computing such metrics, and as we 
show in sequel, we obtain novel properties of the information metric for a wide 
variety of systems, consistent with the physics near the fixed points of these. It 
should be kept in mind that if, in general, the coordinates on the coupling con¬ 
stant space are (x,y), while linearizing about a non-trivial fixed point (x*,y*), it 
is more natural to use as coordinates (Sx, 5y) — (x — x*, y — y*). This should be 
understood by the context. 

Viewed in this perspective, this method bypasses the standard requirement of 
the knowledge of the equation of state (for classical systems) or the many body 
ground state (for quantum systems). Of course this is assuming that the system 
of equations generated from Eq.(j4]) are solvable, but we will show in sequel that 
this is true in a variety of examples. An objection could be that even if Eq. (l4l) 
does yield a solution up to some level in RG, this might not be the case when 
higher order terms are introduced. This is certainly a drawback, which we will 
address towards the end of this paper. For most of this paper, we will concentrate 
on cases where Eq. fjTjl admits an analytic solution. We will see in sequel that we 
are able to capture a large class of models within this simple minded approach. 

As mentioned earlier, important ingredients in any geometric setup are scalars, 
which are invariant under coordinate transformations and indicate global proper¬ 
ties of a curved manifold. In contrast, tensor components such as metric compo¬ 
nents will change under a coordinate transformation, and their scaling relations 
will not in general be coordinate invariant. Ruppeiner conjectured [I] that near 
criticality, R ~ The arguments of Ruppeiner are based on the notion of rela¬ 
tive flatness in a curved space (see Eq.(4.76) of pQ and the arguments preceeding 
this equation). From an RG perspective, we will prove that this result is exact 
up to linear order, but that there are important subtleties when one includes a 
class of higher order terms. Further, the infinitesimal geodesic distance d\ along 
a curve, defined (see Eq.flT])) as d\ = \fg a bdx a dx b is an interesting quantity and 
is known to be related to the concept of a statistical distance. We study this 
object, and show that the geodesic distance is related to the length scale of the 
problem. 

It is well known that that geodesics converge (or diverge) at singularities of 
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a given manifold [21 j. From a celebrated equation due to Raychaudhuri (see 
e.g [22]), the convergence of geodesics can be quantified in this case by a scalar 
expansion parameter (called 0 in sequel). In two dimensional Euclidean cases 
which are of our interest here, the expansion parameter is also an universal in¬ 
dicator of phase transitions, as is the Ricci scalar. There are algebraic relations 
between the three scalar quantities mentioned above, defined on a two dimen¬ 
sional manifold. These give rise to the so called geometric exponents [9]. While 
in that paper, these exponents were calculated in solvable systems, we show here 
that they emerge from the perspective of RG, even beyond linear order. 

This paper is organized as follows. In section 2, we elaborate on linearized 
RG flow equations, generalizing the work of m , and providing a number of new 
results. In section 3, we study non-linear RG flow equations and their geometric 
significance, including logarithmic corrections. Section 4 ends this paper with 
our conclusions and directions for future research. 

Before we embark on our analysis, a word about the notations and conventions 
used in the paper is in order. We consider a variety of examples, and using 
different symbols for the variables will unnecessarily clutter the notation. We will 
proceed with the understanding that the notations used in a particular subsection 
of this paper find usage only in that subsection and are not to be related to 
the other subsections of the paper. Also, the examples used in this paper are 
standard and can be found in textbooks [23], [24]. We will refrain from a detailed 
discussion of the models themselves, which will make the paper unnecessarily 
lengthy, and instead refer the reader to these excellent texts for more details. 
Another important issue should be kept in mind. In order to have a valid notion 
of geometry, the line element d\ 2 = g a bdx a dx b should be positive dehnite. This 
means that along with the diagonal element, the determinant of the metric tensor 
should be positive (see, e.g the discussion around Eq.(3.20) of [1]). In all the 
examples considered in this paper, we have checked that this condition is satisfied. 
We will not mention this in sequel. 


2 Linearized RG and information geometry 

In this section, we demonstrate the construction of the information theoretic 
metric near criticality for linearized RG flows. These may arise in any statistical 
system, when one linearizes the RG equations near a critical point. Here, we will 
be concerned with two parameter examples, and comments on the generalization 
to higher dimensional manifolds will be given towards the end of this section. 

2.1 Case I 

In this subsection, we first recast the results of HH in a form that will be useful 
for us to make some general statements regarding two parameter information 
theoretic models in linear RG, and then go on to study geodesics for these models. 
We start with a theory with two coupling constants x and y, and assume that 
near a generic fixed point ( x*,y *), the linearized RG equations can be written in 
terms of the eigenvalues a and b (these should not be confused with the coordinate 
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labels) as 

x = ax, y = by, (5) 

where the overdot indicates a derivative with respect to a logarithmic length 
scale, l = ln(L). As pointed out in the introduction, appropriate coordinates on 
the PM in this case are (Sx, 5y) = (x — x*, y — y*), which we will still call (x, y) 
by a slight abuse of notation. Let us denote the metric components on the PM 
by g xx , g yy and g xy [25]. 

As mentioned in the introduction, for this example the homothetic vector field 
has components K a = (ax, by) [[26]. We take this as an input and insert it in 
Eq. ([J]) (or, alternatively use K ay> + A^ ;a = Dg a b after lowering the indices of K a ). 
Writing out the components of Eq. ([4]) then leads to differential equations for the 
metric tensor. For the chosen off-diagonal form of the metric, Eq.flU) gives rise to 
the following three equations (as before, a comma denotes an ordinary derivative 
with respect to the variable that follows it) : 

gxx,x A byg xx y A (2 a D)g xx 0 , 
axg y y : x A byg yy y -1- (26 D)g yy 0 , 

(^xgx y ,x A byg xy y A (a A b D)g xy 0 . (6) 

This reveals that the metric components are generalized homogeneous functions 
near criticality mi and immediately reproduces the well known static scaling 
relations. This is true for any linearized set of RG equations. Say the variable x 
drives the phase transition. Then we can write the general solution for Eq. ()6|) as 

gxx = x^~ 2 Gi [yx~^j , g yy = x 3 ~^ L Q 2 [yx~^j , g xy = x~ B±i ^G 3 [yx~^j . (7) 

Here, Gi, i = 1,2,3 are functions of a single variable yx~ and reminiscent of 
Widom scaling of the free energy near criticality [23], |24]. However, the functions 
Gi are not the same as scaling functions that apper in the free energy, since the 
metric is more naturally interpreted as the second derivative of the free energy 
in the conventional picture of information geometry. It will be assumed that the 
functions G% are analytic and equal a constant value near criticality, y = 0. We 
will advocate two arguments to justify this, and go on to check these with the 
known example of the 1-D Ising model. 

First, note that in the classical notion of information geometry, the metric 
components are related to the classical response functions. For example, in a 
magnetic system, if x is identified with the reduced temperature t = (T — T c )/T c , 
with T being the temperature and T c its critical value, and y is identified with 
the reduced magnetic field H/T c , then the metric components g xx , g yy and g xy 
are related to the specific heat, the magnetic susceptibility and the derivative of 
the magnetization, respectively. Eq. (EJ) then indicates that these have the correct 
critical exponents if the functions Gi are analytic at y = 0, and equal a constant 
of order unity. 

This can also be seen by noting that Eq. (l6]l translates into the fact that the 
metric components are generalized homogeneous functions up to linear order so 
that standard scaling arguments can be applied (see Eq.(5.5) of [II]). This is 
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indicative of the fact that the functions Q h can be taken to be analytic, and 
constants of order unity close to criticality. There is a small subtlety here. If 
we take all functions Q t to equal the same constant near criticality, the metric of 
Eq.(]7]) become singular. Hence, this should be avoided, and Qi , 2 = 1,2,3 have to 
be taken to equal different constant numbers of order one. These multiplicative 
constants can at most affect our results for the metric, the scalar curvature and 
the expansion parameter by some overall constants, and will not affect our scaling 
analysis. Keeping these explicitly in the computations will clutter the notation, 
and without loss of generality, we will take two of these to equal unity, and set 
the other one to 2. This is just a particular choice and any other choice would 
only affect the results by an overall numerical constant. With this choice, we are 
also able to consistently satisfy the positivity constraint on the line element, as 
can be checked. 

Secondly, assuming that the scaling function is analytic, and equals a constant 
of order unity near criticality, we obtain the scalar curvature in terms of the 
driving parameter x up to an overall constant as 

o h 

R = —(D-2b)x~^ (8) 

a 2 

This shows that the scalar curvature blows up if x is relevant, and goes to zero 
of x is irrelevant. If we assume that R ~ £ D where £ is the correlation length, 
then Eq. (JH]) implies that £ ~ x~^, i.e correctly reproduces the correlation length 
exponent for classical phase transitions where x is identified with the reduced 
temperature. For quantum phase transitions, we note that if a: is a relevant 
variable, then a perturbation in this direction produces a gap in the spectrum 
that in turn indicates that the correlation length scales as x ~+ This is again 
consistent with R ~ £ D (and also justifies the assumption of R ~ £ D ). We 
thus see that assuming that the functions Qi are constants of order unity near 
criticality produces a consistent geometric picture up to linear order in RG. This 
will be assumed in what follows. We will not explicitly indicate these functions 
in sequel. 

ft is instructive to validate our analysis thus far by comparing it to a known 
example. We choose the standard example of the classical 1-D Ising model in a 
magnetic held, with the Hamiltonian given by 

N N 

H = -j'£S i S j »-h'£ S > < 9 > 

3 =1 3 =1 


Information geometry for this model was worked out in [[27] in the limit of large 
N, and we quote their result for the metric. First we define the variables x = J/T 
and y = h/T (where we set the Boltzmann’s constant to unity). Further, writing 
t = e~ Ax and near the critical point substituting t = e and h = 5 , the metric 
components (Eq.(4.17) of [27]) read, after some algebra, 


Qxx 


e 2 1 + 2(he 2 ) 2 

4 (l + (5e-i) 2 ) i 


9xy 


e -1 8 1 _ _i 1 

~n~ ~ T TT ’ 9yy ~ e 2 + 7 

62 (l + (he4) 2 ) 2 (l + (<5e^) 2 ) 

( 10 ) 
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Thus the metric is similar to the one in Eq. flTj) (with x = e and y = 8) upon 
identification a = 2, b = 1, with D — 1. The functions defined in that equation 
read 


Gi 

(&T*) 


1 l + 2(5e-i) 2 

^(l + (£e ^) 2 ) 5 

1 

(l + (<5e^) 2 ) 5 


£2 (5e-s) 


2 


5e 2 




( 11 ) 


which are analytic near criticality, as expected, if we assume e and <5 to be of the 
same order. Note that these are different functions which evaluate to different 
numerical values near criticality, as alluded to before. Also, from our discussion 
it follows that the RG flow equations here are governed by t = 2t, h = h, in 
agreement with Eq.(32) and (33) of [18] . The analysis of the scalar curvature 
and geodesics for the one dimensional Ising model has been done in [9j, to which 
we refer the reader for more details. 

The metric of Eq. (]7|) is to be used when x is the driving parameter in the 
phase transition. An equivalent form of writing the solutions of Eq.flH]) is 


lj — za / _ a \ lj — zo / _ a \ _ o.-|-b — D / _ t£ \ 

= V * (-''!/ *), 9m=V 1 V. [xy 1 ) , 9:„1 = ;/ * Vj (.ri/ ») . 


( 12 ) 

Here J-), % — 1, 2, 3 are arbitrary functions of the variable xy *> and are assumed 
to approach a constant value at x = 0. The metic of Eq. ffl2il should be used 
when y drives the phase transition. This metric has a scalar curvature given by 


R=^D-2a)y-% (13) 

Comparing Eqs.ffSj) and (TT51) . we see that the divergence of the scalar curvature 
is controlled by the coefficient of the driving parameter in the RG, which is 
expected. As before, the functional forms of T t (xy~%) may be different for 
i — 1,2, 3. However the only assumption here is that all these are of order unity 
and that the leading behaviour of the metric near criticality is controlled by the 
exponents of y. 

Eq.IjSj) is applicable to any model of linearized RG and shows that if the 
parameter x is relevant, the divergence of the scalar curvature is controlled by 
the relevant eigenvalue. If we identify the scalar curvature with the correlation 
volume (up to a possible arbitrary constant), it is seen that the correct correlation 
length exponent is recovered. 

An interesting quantity that we will now focus our attention on is a set of 
geodesics (called a geodesic congruence) on the PM. Geodesics are analogues of 
straight lines in curved spaces, and these are paths that minimize distances be¬ 
tween points on a curved manifold. For curved information theoretic manifolds 
that we describe here, geodesics provide a further characterization of classical and 
quantum phase transitions, as shown in [28] . Namely, one considers a geodesic 








congruence on the PM, and it can be shown that near a critical point, the con¬ 
gruence converges (or diverges). 

Let us make this statement more precise. If our PM is defined by the coor¬ 
dinates (i.e coupling constants) x a , then geodesic paths on the manifold satisfy 
the equation ( x a )" + T^ c (x b )'(x c )' = 0. Here, = \g ad (, g db , c + g dcfi - g bc4 ) are 
the Christoffel connections defined from the metric, and the prime denotes a 
derivative with respect to an affine parameter A along the geodesic, which is con¬ 
ventionally taken to be the square root of the line element, i.e dX 2 = g ab dx a dx b . 
For such an affinely parametrized geodesic, the geodesic equations can be ob¬ 
tained from a variational principle from the Lagrangian £ = | (g ab (x a ) ( x b ) ). 
This will be useful for us later. 

If we denote as the normalized tangent vectors u a = ( x a ) , curvature effects 
on geodesics (near criticality) are measured by the tensor B a h = V b u a . Recall 
that the covariant derivative on a generic vector V a is defined by the action 
V a V b = d a V b + Y b a( V c . Then it can be shown that 0 = B a a , called the expansion 
scalar, gives an effective measure of the convergence or divergence of a geodesic 
congruence. At critical points on the PM, i.e at phase transitions, 0 diverges. 0 
being a scalar quantity, this is a coordinate independent characterization of phase 
transitions. To compute 0, we require a solution for the vectors u a . In general 
this might be difficult to obtain, but when the metric is independent of one of 
the coordinates (as will always be the case here), such solutions can be found 
analytically if u a is normalized, i.e u a u a = 1. These analytic solutions constitute 
our geodesic congruence. We refer the reader to [9] for more details. 

To illustrate the procedure, we first recall the two dimensional metric of 
Eq.flTJ), where we will include a multiplicative factor of k\ in the g yy compo¬ 
nent. Now denote the tangent vectors to geodesic trajectories by the vector 
u a = (x'{X),y'(X)), where a prime denotes a derivative with respect to an affine 
parameter A. Normalization of u a imposes the condition 


/ \ \ —2a — 2b+D 

x[X) 


x\X) 2 x(\) « + 2x'(X)y' (A)x(A) “ + kix(X) 2 y' (A) 2 


1 (14) 


The left hand side of Eq. (fl4ji is in fact proportional to the Lagrangian alluded 
to before. Noting that this is independent of the coordinate y (as is the metric), 
the Euler-Lagrange equation for y( A) imposes a further constraint dC/dy' = fc 2 , 
where P 2 is an arbitrary constant. Then, if this constraint is solved in conjunction 
with the normalization condition, we obtain 


x'(X) = Jx{\) 2 ^2x(A)« — k^x^X) 2 a j, 

y'{ A) = ^ | k 2 x(X)—^— — x(X)~^~ 1 x(X) 2 ^ 2 ^ ^2x(A)? — fc|x(A )^j j (15) 


The first of these equations can be solved to obtain an expression for the geodesic 
distance in terms of the Gauss Hypergeometric function, 


A 


\/2 D_ 

——ax 2a 2 Pi 


1 D D 1 2 26-d\ 

2 ’ 46 — 2D ’ 46 — 2D + lj 2‘^ ° ) 


- hi 


(16) 
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where k 3 is another arbitrary constant. I 11 order to obtain a real value of A which 
is physically reasonable, we require b > -D/2, and a further constraint on the 
constant k 2 . For small values of x, this can be seen to restrict the Hypergeometric 
function to values close to unity. This last fact indicates that it is reasonable to 
set k 2 = 0, without loss of generality. Indeed, Eq. ffThl) simplifies in this limit 
and we obtain as a solution, A ~ x ^ — k 3 . Now note that we are interested in 
geodesics that reach very close to the critical point. It is natural to measure A 
from the critical point, so that we will require A —» 0 as x —> 0. This indicates 
that the constant k% can be set to zero as well [29]. 

From Eq. (ll5p . we now obtain the following solutions for x and y as a function 
of A : 


x(A) = 2 


a 

D 



y( A) 


2- 


b+P 

d a 


' D\\ 


2b 

D 


b 



(17) 


where we have imposed the condition x{\ = 0 ) = 0 , i.e the affine parameter 
is measured from criticality. These equations can be now inverted to obtain an 
analytic expression for the affine parameter, namely, A ~ x^, apart from constant 
factors. Also, using Eq. dlSP and the metric of Eq. (J7]), we obtain by some simple 
manipulations, 

(-D — 2 b)x~^ 

" = ~ VTa 



Using the solution for the affine parameter, we obtain as x —> 0 from Eqs. (fT3l) 
and (HEP, 

-R ~ A 2 , 0 ~ A 1 (19) 


The same conclusion can be reached for the metric of Eq. (fT2lh as can be easily 
checked. Eq. dlfJD can be thus understood as an universal indicator of phase 
transitions for any system with linearized RG flow. The exponents appearing in 
this equation were dubbed as geometric critical exponents in [9]. In that paper, 
the analysis was conducted by exploiting the behavior of the information metric 
close to criticality for exactly solvable systems. Here we have given proof that 
the relations hold for any arbitrary two parameter system, at least up to linear 
order in RG. 

We record a couple of observations before we move on. From Eq. (1T31) and 
(in, note that the scalar curvature and the expansion parameter diverge if the 
operator x (or y) are relevant. If these are irrelevant, i.e a or b is negative so that 
x or y are stable directions, then these go to zero in the limit that the coupling 
constants go to zero. In that case, A calculated from Eq. (fm (after setting k 2 and 
to zero) approach infinity. This is a typical feature of information geometry 
that we will come across later also [30]. Note that the relations of Eq. (TT9l) remain 
valid irrespective of whether the operators are relevant or irrelevant. 

Also note that the geodesic distance A can be related to the length scale of the 
problem as follows. Using the fact that under an RG transformation, £ = £oe*, 
we End that R ~ £ D translates into l = —(2/D>)logA. This is a generic feature 
for all linearized cases. 
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2.2 Case II 

Our next example is that of an RG flow of the form 

x = a\x + a 2 y, y = hy. (20) 

This form of the RG equations occur in perturbation theory, for a linearized one- 
loop approximation in Landau-Ginzburg models. In this case, one can obtain a 
metric using the original variables, as we illustrate in a moment. The important 
point is that geometric methods can be applied to set of redefined coordinates, 
consistent with the eigendirections of the RG flow equations. For example, in 
this case, if we define a new variable 0 = x + a 2 y/ (cq — 61 ), the RG flow equations 
reduce to z = cqz, y = byy. The results of our previous analysis can now be 

readily applied in this new set of coordinates. In particular, for y = 0, we obtain 

the components of the information metric as 

D o D — 2bi ai+bi—D 

9zz = Z ai , 9yy = Z »i , g yz = Z “1 (21) 

and an entirely similar analysis holds for z = 0. In both cases it can be seen 

that our previous results R ~ A ~ 2 and 0 ~ A ^ 1 hold. As before, the scalar 
curvature can be computed and the expressions are similar to the ones in Eq. () 8 l) 
and m and with the identification R ~ £ D shows that while oq is the critical 
exponent in the direction y = 0, b\ is the one in the direction z — 0. We have 
thus constructed the information metric for the Landau-Ginzburg model at one 
loop, Eq. (j 2 T]h solely by using the RG flow equations. 

We should mention here that using the original set of equations (Eq. (|20l) ) 
also, it is possible to compute the metric tensor. However, this has a complicated 
structure and does not reveal any meaningful physics. Our scalar relations are 
however expected to hold here as well. It should thus be kept in mind that to 
interpret the various quantities associated with information geometry, one needs 
to correctly choose coordinates. Once this is done, analysis of the metric be¬ 
comes meaningful, and the correlation length exponent comes out correctly with 
the identification R ~ with D being the spatial dimension of the system. To 
summarize, for any set of linearized RG flow equations, the information metric 
can be written down simply from the scaling dimension of the operators. An ap¬ 
propriate choice of coordinates then predicts the correct exponents of the system. 

2.3 Case III 

Before we close this section, we will comment on a situation in which a system has 
a critical line, for example a gapless line in the parameter space for quantum phase 
transitions. This is exemplified by the one-dimensional anisotropic Heisenberg 
spin 1/2 chain. This model was considered in [32] . where the gapless line was 
interpreted as a spin flip transition. The model Hamiltonian is given by 

H = E K 1 + 7)s;s“ +1 + (1 - 7 )S’S’ +1 + A S-S* +1 - hS a (22) 


11 




where S l , i = 1,2,3 are spin operators, 7 is an anisotropy parameter, A is 
the coupling in the ^-direction and h is a magnetic field along the x direction. 
As shown in [32] . bosonization techniques yield the following perturbative RG 
equations in terms of h, 7 and b (the latter being the coefficient of an operator 
which arises in an operator product expansion) are 

h = aih — a 2 ^h — a 3 bh, 7 = 617 — b 2 h 2 , b = c\b + c 2 h 2 , (23) 

where the coefficients are determined by the scaling dimension of the correspond¬ 
ing operators and read 

ai = 2 - K- 1/(410, 02 = 2 -1/K, a 3 = 2-4/1 
b 1 = (2 — 1/ K), b 2 = (2K — 1/(2K)) = c 2 , Cl = (2 - 4 K) (24) 

Here, K is related to A (Eq.(3) of (32] ) and takes values 1/2 < K < 00 (for 
details, see [32]). Now suppose we are at a fixed point h = h*, and look at 
information geometry in the 7 — 6 plane, where the RG equations are linear. We 
will not go into the details here, but simply state the result that the information 
metric at a fixed point of h can be obtained to be 

n K l-K(lK-3) K(3—8K) 

grr = z " , 9"/b — z , g bb = z i-** (25) 

where z = (2y(l — 7 ^) — 2h* 2 (K — -j^)), and a coordinate defined by y — 6(2 — 
4 K) + h* 2 (2K — 1/(2 K)) is set to zero. The scalar curvature and the expansion 
parameter diverge as 

2 K K 

Rr^z 2K ~1, 0 ~ Z 2X -! (26) 

and it can be checked that R ~ A^ 2 and 0 ~ A -1 , with A being the geodesic 
length, as expected. For different values of K within its specified range, the 
scalar curvature and the expansion parameter diverges at 7* = (2 K + l)h* 2 /2. 
An entirely similar analysis holds for z = 0 when the information metric is 
determined by y. This results in b* = y*/(2 K). These values of 7* and b* 
determine a fixed line in the 7 — b plane and also determine the value of h* which 
is entirely consistent with [32], proving the validity of Eq. ([25]) . 

As an aside, we point out that a similar RG flow equation as in Eq. ([23]) was 
obtained in [[33] for a model of two weakly coupled Luttinger chains, where the 
second and third terms of the first of Eq. ([23]) were absent. In this case, a simple 
transformation of variables £1 = 7 — b 2 h 2 /( 2 cq — b\) and z 2 = b — c 2 h 2 / ( 2 ai — < 7 ) 
renders the RG equations linear, i.e h = aih, Z\ = biZ\, z 2 = C\Z 2 . These 
equations can be used to define the homothetic vector field in three dimensions, 
where an entirely similar analysis as our two dimensional examples so far can be 
done. We will however not present the details here : higher dimensional examples 
will be treated elsewhere. 

Finally we point out that setting A = 0 in Eq. (j22|) reduces it to the transverse 
XY model considered in the context of information geometry in [5]. This model is 
exactly solvable, so as a curiosity we check what our RG method predicts as the 
information metric. A bosonization procedure near 7 = 0, h = 0 here yields the 
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RG equations 7 = 7 , h = h. Along with Eq. flT]) and the fact that for linearized 
RG flows, the information metric can be taken to be diagonal [25], we recover 
g 77 ~ 7 _1 , gth ~ 7 -1 , in agreement to the metric derived in Eq.(7) of [5], near 
7 = 0. 

3 Nonlinear Examples and Logarithmic Correc¬ 
tions 

In this section, we will focus on cases where the RG flow equations include non¬ 
linear higher order terms. As in the previous section, an appropriate combination 
of variables will be seen to render the solutions tractable. 

3.1 Case IV 

Our Erst example is given by the flow equations 

x = aix + a 2 x 2 , y = biy (27) 

These flows arise, for example, in the critical dynamics of a time independent 
random held interaction introduced in a Ising spin or quantum rotor model |31] 
in the e expansion, above the lower critical dimension d c = 2. The Hamiltonian 
for the model is 

H = (28) 

(*?> * 


where T is the strength of a transverse held, and hi are site dependent magnetic 
helds, with the random helds correlated in one direction. RG hows in this model 
were studied in [3T] and the equations are of the form of Eq. (l27lh with the iden¬ 
tification x — h/J and y = T/Jq where T is the temperature, h is a measure of 
the randomness of the random held hi defined via its distribution (Eq.(2) of [ 8 l]). 
Jq is the interaction in the direction in which the helds are correlated. Also, for 
this model, aq = —e /2 and fq = — (1 + e), with d = 2 + e. In general, obtaining 
the information metric for such a system is not possible exactly, but can be done 
using the symmetry of the RG how equations, as we show below. 

To keep the discussion general, we proceed with arbitrary cq, a 2 and fq. First, 
note that a transformation of variables z = x/{a\ + a 2 x ) renders the equations 
linear in the variables 2 and y, i.e i = (i\Z and y = b\y. Now we can construct 
the information metric using the methods described in the previous subsection. 
Since both x and y are irrelevant up to linear order, i.e are stable directions, 
following our previous discussion the scalar curvature does not blow up here. 
This is the trivial hxed point. In order to obtain the information metric near 
the non-trivial hxed point, it is more convenient to work in terms of the original 
variables, and we illustrate the results below. The differential equations for the 
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metric components are obtained in the variables x and y as 


2(cii “I - 2o 2 x^g xx T b\yg xxy T x(cl\ T d 2 x^g xxx Dg xx 0, 

(«1 T bi 4- 2d2X^g xy T biyy X y,y T ic(oi 4- x)y X y,x Dg xy 0, 

2b 1 g yy + b x yg yy , y + x(a 1 + a 2 x)g yy)X - Dg yy = 0. (29) 


These have the solutions given by 


9 XX 
9yy 


D o _ D o aid-bj—-D a l~ b l+^ 

x°i “ (a 2 a; + ai) a i , g xy = x °i (a 2 x + cp) a i , 




“I - 


D-26-1 

a l 


(30) 


The geodesic equations can be obtained in the same way as outlined in the pre¬ 
vious section, and denoting the tangent vector as (x'(A), y'(X)), we obtain as a 
normalized solution 


x'(X) = V2x(ai+d2x) 


x 


d\ + & 2 X 



if X 
y/2 \&l + CI 2 X 


D 

2a^ 


a l 


(31) 


The above equation then yields 

\pl / X 

D \a 1 + a 2 xj 


(32) 


Noting that the scalar curvature and the expansion parameter are given by 


R 


2 h (D 


26i) 


x 


CL 2 X + d\ 


D 

a ! , 0 


D ~2bi f x A __g_ 
y/2 \a 2 x + aij 


(33) 


we obtain R ~ A~ 2 , 0 ~ A -1 , as expected. Note that these are true in D = 
2 + e dimensions. Since d\ = —e /2 is negative, the curvature and the expansion 
parameters diverge at the non-trivial fixed point x* = —cq/cp and t* = 0. There 
is no divergence of these quantities at the trivial fixed point x* = t* = 0 . 


3.2 Case V 

Next, we come to the case where the RG equations are taken to be 

x = a\x + a 2 x 2 , y = biy + b 2 xy (34) 

These is the standard textbook form for RG equations of Ising like models near 
four dimensions, i.e in d — 4 — e, if we identify d\ = e, a 2 = —72, b\ = 2, b 2 = —24. 
For O(N) vector models, the coefficients a 2 and b 2 are given by —8(n + 8) and 
—8(n + 2 ) respectively [24] . 

First note that upon making the coordinate transformations Z\ = xj(di+a 2 x), 
b 2 

z 2 = y{a 1 + d 2 x) a 2 , these equations reduce to £\ = a±zi, £2 = b\Z\. This clearly 
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defines the Gaussian fixed point, with the critical exponent along the y direction 
being b\ = 2. Next, if we linearize the RG equations near the Wilson-Fisher fixed 
point x* = —ai/a 2 = e/72, we get back the thermal exponent 2 — e/3, as follows 
from our discussion of the previous section. 

ft is interesting to see what the information metric reveals when applied to 
non-linear set of equations, Eq. (l34lh Here, it is convenient to first define a new 
variable z = y(a 2 + ai/x) bl ^ ai , so that we have an equivalent set of equations 

x = a\x + a 2 x 2 , z = b 2 xz (35) 


Note that since b\ and e are positive, the coordinate transformation mentioned 
above becomes ill defined near x —> 0, i.e the Gaussian fixed point. However, 
this is not the case near the Wilson-Fisher fixed point. We will need this fact 
later. We take the components of the homothetic vector held, generating the 
scale transformation near criticality as (x,z). Then, the equations determining 
the metric on the parameter manifold (with parameters (x,z)) are given by 

x ( 2b 2 g zz + b2zg ZZ)Z + (cq + a 2 x)g ZZiX ) - Dg zz = 0 , 

(ai + x(2 a 2 + b 2 ) g xz + b 2 zg zz + x(b 2 yg XZiZ + («q + a 2 x)g XZtX ) - Dg xz = 0 , 

2(tq T 2a 2 x)g xx T 2b 2 zg xz T x ( b 2 zg xx ,z T (cq T O j 2 x)g XXiX ) Dg xx 0. (36) 


The metric can be determined from the above equations by first solving for g zz 
and hence g xz and g xx . The expressions are lengthy, but simplify in the limit 
z = 0 (equivalently y = 0 ) and read 

D o _ D_ _ 2 _D__^ &2 a\ + D 

9xx = x ai (a 2 a; + ai) »i , g xz = x a i (a 2 x + “2 «i , 

JL 2b 2 D 

9zz = x a i ( a 2 x + ai) “2 °i. (37) 


Having obtained the metric, we now focus on geodesics on the PM. As usual, 
we start with the tangent vector (a:'(A),P(A)). The fact that the metric is inde¬ 
pendent of z leads to 


x'(A) = V2.X 1 2o i (ai + a 2 x ) 1+2a i , z\ A) 



D 

2 


_£-^_ _J_ 

(eq + a 2 x ) 2a i 


b 2 

a 2 


(38) 


The hrst of the above equations can then be solved to give 


y/2 f x A 2 “i 
D \ai + a 2 xj 


(39) 


which was also the expression obtained for the quantum rotor model studied in 
subsection 3.2. The scalar curvature and the expansion parameter are given by 


R 


J-J LJ 

—2b 2 x °7 (a 2 x + ai) “7 (—2 a 2 x — 2a\ + 2 b- 2 x — D ) 

2 =(d- %) (— 

y[2 V°i +a 2 x) 
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0 


(40) 





Now note that in this case, a\ = e is positive. Hence, the scalar curvature 
and the expansion parameter seems to diverge at the Gaussian hxed point x —>■ 
0. However, as mentioned before, the coordinate transformation used to derive 
Eq. (l35l) is not trustable here. Let us thus focus on the Wilson-Fisher hxed point, 
where from Eq. (l39lh we also see that A —> oo. This may look at variance with the 
linearized result (where the geodesic distance goes to zero near a non-trivial hxed 
point), but that is not so. To see this note that in our non-linear analysis, we 
have not linearized about the Wilson-Fisher point. The geodesic distance here is 
still measured from x = y = 0 , so our analysis simply rehects the fact that the 
Wilson-Fisher hxed point is an infra-red hxed point. By combining Eqs. (13911 and 
(flOl) . we recover the relations R ~ A -2 and 0 ~ A -1 at this hxed point. 

3.3 Case VI 

Another interesting situation occurs when one of the variables is marginally ir¬ 
relevant. We will focus on a set of RG how equations of the form 

x = cliX 2 , ij = biy + b 2 xy (41) 

These are the RG how equations for the Ising model in four dimensions (with 
ai = —72, b\ — 2 and b 2 = —24), as follows from our analysis of the previous 
section, by setting e = 0. Clearly, the results of that analysis cannot be used 
here simply by setting the coefficient of the linear term to zero, since the metric 
components of Eq. ()37l) are ill defined in this limit. To perform this analysis, we 
start by defining a new variable z = ye bl ^ aiX \ Since a\ is taken to be negative, 
and hence this transformation is well defined in the limit x — > 0. It follows that 
the RG equations can be written more conveniently as 

x = aix 2 , z = b- 2 xz (42) 

In terms of the variables x and z, the equations for the metric components in this 
case are seen to be 

x (2 b 2 g zz + b 2 zg ZZtZ + a ± xg ZZjX ) - Dg zz = 0 , 

(2ai + b 2 )xg XZ + b- 2 zg ZZ + x ( b 2 zg xz,z “h CL\ %Qxz,x) - Dg XZ 0 , 

Ao,\xg X x — 2 b 2 zg X z x ( b 2 zg XXjZ T ci\xg XXtX ) T Dg xx 0 . 

As before, the hrst of these equations can be solved to obtain g zz , which can in 
turn be used to fold g xz and hence g xx . The solutions for the metric components 
read 

1 _2f2_4_ D_ / &2 

9xx = -oX ai e “ 1 - l ai x a i 

aj V 

2&2 D 

9zz = hx “1 e “ 1 " 

Here k\ is a constant that we will set to 2 following our previous discussions. We 
hrst record the expression for the scalar curvature calculated from the metric of 


b 2 xz) , g x 


1 2 (a 1 +b 2 ) D / _^2_ 

—x a i e a i x (ax a i — b 2 xz 

a\ \ 


( 44 ) 
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Eq.flTTl) in the limit z = 0, 


D 

„ 2b 2 xe a ^ x (2a x x - 2b 2 x + -D) ,, rX 

R =-j- (45) 

a i 

Clearly, with negative a\, R rapidly goes to zero for small values of x. Now let 
us consider geodesics on the manifold defined by the metric of Eq. (ITT]) . For the 
four dimensional Ising model, a\ = —72 and b 2 = — 24, and for small x and z, we 
can ignore the xz pieces in Eq. flTTl) . which renders the said metric independent 
of z. As before, we consider a tangent vector to a geodesic, which we denote by 
(x'(A), z'(A)). Here, A is an affine parameter along the geodesic. Normalization 
of the tangent vector, along with the fact that the metric does not depend on z 
implies that 


x'(\) = V2x(X) 2 e 2a ^ w , z'{ A) 



&2 D 

a-^ g 2a^ x(X) 


Solving the first of these equations, we obtain 


x(A) 


D 

2a > lo s (vfir) 


(46) 


(47) 


If is negative, as happens for the four dimensional Landau-Ginzburg model, 
we see that at x = 0, A has to go to infinity, although logarithmically [35]. After 
some algebra we obtain here, 


R = 


2 b 2 


( a i '°g ~ 



In the limit A —» oo, we finally obtain 


R ~ 


1 

A 2 logA’ 


0 r*j 


1 

A 


(48) 


(49) 


The first of these equations indicate that the Ricci scalar picks up logarithmic 
corrections to geometric scaling relations in four dimensions. 


4 Conclusions 

In this paper, we have provided evidence that scale invariance in the vicinity of 
a critical point can provide valuable information on the metric of the parameter 
manifold in classical and quantum phase transitions, in a unified fashion. In 
particular, this method can be applied to systems that are not exactly solvable 
to read off the scaling behavior of the metric (and hence related quantities like 
the fidelity susceptibility) near criticality. Our method complements the work 
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of 0) to determine scaling patterns for information geometric quantities. While 
the work of HU utilises the RG equations up to first order to read off the scaling 
of the metric in classical phase transitions, our method explicitly solves for the 
metric components in a variety of non-linear examples. 

While most of the paper deals with classical phase transitions, we have given 
one non-trivial example of a quantum phase transition that can be studied in 
this framework. This extends the study of information geometry to novel settings, 
which, to the best of our knowledge, have not appeared in the literature. We have 
seen here that the relations R ~ A -2 and 0 ~ A -1 are universal, except for the 
four dimensional Ising model, where the former relation picked up logarithmic 
corrections. This strengthens the claim made in |9j about universal geometric 
critical exponents. 

We should mention here that in two dimensions, the scalar curvature R and 
the expansion parameter 0 satisfy the Raychaudhuri equation 0 2 + 0' + R/2 = 0, 
where a prime denotes a derivative with respect to the affine parameter A [34j . In 
principle, given R (0), this equation can be used to determine 0 ( R ). However, 
it is not always possible to solve this equation analytically. In all the examples 
considered in this paper, as a cross check on our results, we have verified that 
the Raychaudhuri equations are indeed satisfied. 

In this paper, we have considered a class of examples where the information 
metric was obtained from the set of RG flow equations. Clearly, one might argue 
that this may not be the case for more generic examples. Consider for example 
an RG equation of the form x = a±x + a 2 y, y = b±y + b 2 x. In this case, the 
homothetic vectors do not have an analytic solution as can be checked. This is 
a caveat in our analysis. A further criticism might be that higher order terms 
than those considered here, are difficult to take care of. We note however that 
in general, such terms might be solved iteratively, i.e one can use perturbation 
theory to solve the differential equations for the components of the metric tensor. 
This is substantially more complicated that the analysis presented here, and a 
full study of the same is left for the future. 

It will be interesting to extend the present analysis to cases where the param¬ 
eter manifold has dimensionality higher than two. One example was commented 
upon in this work, but a broader analysis might reveal interesting facts about the 
geometry of the renormalization group, as higher dimensional PMs offer more 
structure, and in particular, more scalar invariants. What these scalars mean in 
the context of RG will be an interesting issue for future investigation. It might 
also be interesting to consider the role of time in information theory [36], in the 
context of the models considered here. This issue is currently under investigation. 
Finally, it might be useful to investigate information geometry in the context of 
Kosterlitz-Thouless type phase transitions in the two dimensional XY model. 
Preliminary analysis indicates that here the scalar curvature of the information 
metric diverges exponentially, in lines with the behavior of the correlation length. 
However, this case requires further understanding. 
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